Discrete Thermodynamic Bethe Ansatz 



Michel Bergere + , Ken-Ichiro Imura ++ and Stephane Ouvry ++ 

+ Spht, CEA-Saclay 
++ LPTMS, Universite Paris-Sud 
bergere@spht.saclay.cea.fr; imura@ipno.in2p3.fr; ouvry@ipno.in2p3.fr 

(February 1, 2008) 

Abstract 



We propose discrete TBA equations for models with discrete spectrum. 
We illustrate our construction on the Calogero-Moser model and determine 
the discrete 2-body TBA function which yields the exact iV-body Calogero- 
Moser thermodynamics. We apply this algorithm to the Lieb-Liniger model 
in a harmonic well, a model which is relevant for the microscopic description 
of harmonically trapped Bose-Einstein condensates in one dimension. We find 
that the discrete TBA reproduces correctly the iV-body groundstate energy 
of the Lieb-Liniger model in a harmonic well at first order in perturbation 
theory, but corrections do appear at second order. 
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I. INTRODUCTION 



It has been known for a long time that the spectrum of the Calogero model [ffl], defined 
here as particles on an infinite line interacting via a(a — Xj) 2 2-body interactions, can 

be found by the Bethe Ansatz (BA) which assumes periodic boundary conditions, and that 
its thermodynamics can be obtained, in the thermodynamic limit, by the Thermodynamic 
Bethe Ansatz (TBA) ||. It is however not necessary to rely on the BA (or the TBA) to 
get the Calogero spectrum (or its thermodynamics). Indeed, the Calogero model is exactly 
solvable 

- either by confining the particles in a harmonic well of frequency u: the Calogero- Moser 
model Q with discretized harmonic well quantum numbers and energies, and Hamiltonian 

Hn = ~\ E ^ + a(a - 1) £ + \u 2 x 2 (1) 

- or by confining the particles in a periodic box of length L: the Calogero-Sutherland 
model || with discretized momenta and energies, and Hamiltonian 

Hn = -\ E f~ 2 + <*(<* - l)(y) 2 E V] (2) 

(the 1/ sin 2 [7r(xj — Xj)/L] interactions are nothing but the periodic version of the infinite line 
interactions). It is not a surprise that Bethe ansatz equations yield the Calogero-Sutherland 
spectrum, since they also assume, as stressed above, periodic boundary conditions. 

Both parameters u and L can be considered as long distance regulators, the thermody- 
namic limit, i.e. the infinite line limit, being obtained either by uj —>■ or L — > oo, resulting 
in continuous momenta and energies. 

The Calogero model describes particules with intermediate statistics, which is natural 
due to the topological (statistical) nature of the l/(x, — Xj) 2 interaction in Id. In the ther- 
modynamic limit indeed 0, the Calogero thermodynamics realizes microscopically Haldane 
(Hilbert space counting) statistics 0. Moreover, the Calogero model has been shown to be 
obtained as the vanishing magnetic field limit |§ of the lowest Landau level anyon model 
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(LLL-anyon model) in the regime where the flux tubes carried by the anyons screen the flux 
of the external magnetic field (screening regime). Not surprisingly, the LLL-anyon model 
also realizes microscopically Haldane statistics, the Hilbert space counting argument being 
manifest here via a mean field argument (adding anyons screen the external magnetic field, 
and thus diminish the Landau degeneracy of the total -mean+external- magnetic field): thus 



a clear relation between Haldane |10[ and anyon statistics |11 . 

Starting from the BA spectrum, and following Yang and Yang footsteps |§, one can 
compute a la TBA the thermodynamics of the Calogero model in the thermodynamic limit 
L — > oo. The thermodynamical potential In Z -where Z = J2n=o zN %n is the grand partition 
function- ends up to be those of a Fermi gasQ 

log Z = ^- f°° dk log[l + ze-M% (3) 

but with a 1-body energy e(k) defined in terms of the free continuous 1-body quadratic 
spectrum e (k) = k 2 /2 as 

Pe(h) = Pe (h) - A r dk^{h - k 2 ) log[l + ze~^) (4) 

All J-oo 



In the Calogero case [12 



^{k 1 -k 2 ) = ^{l-a)6(k 1 -k 2 ) (5) 



lr rhere is an equivalent formulation in terms of a free Bose gas, namely 

L f°° ., , 1 



L r°° 

log Z = — / dk lo. 

27T J-oo 



1 - ze~^ 



but with a 1-body energy e(k) defined as 

L r°° i 

f3e(h) = Peo{ki) - — / dk 2 Hh - k 2 ) log 

One has obviously 



1 - ze-P^ 



2tt 

*(fci - k 2 ) - - k 2 ) = —j-Sih - k 2 ) 
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is intimately related to the 2-body scattering angle, and encodes, if one thinks in terms of 
statistics, the statistical exclusion between two quantum states, here with the same momen- 
tum. 

Note that if one denotes y(k) = 1 + ze' 13 ^, which can be regarded, in view of ([!]), as 
the grand partition function at momentum k, then (|j) can be rewitten as 

y{k) - ze-^^yik) 1 - = 1 (6) 

a particular case of Ramanujan equations ||13|| . 

As already said, equations of the type (|3],|6]) were first obtained directly by i) considering 
the exact iV-body Calogero spectrum in a harmonic well 0, or by considering the exact 
iV-body LLL-anyon spectrum in a harmonic well ||, ii) and then taking the thermodynamic 
limit. 

Now we might ask the following question: are the TBA equations (HfD specific to the 
thermodynamic limit with continuous momenta and continuous dressed energies, or can 
they also describe the thermodynamic of the Calogero model in a harmonic well or in a 
periodic box with discretized energies? In other words, can we find a discretized version 
of the function — k 2 ) in ([5|) such that the harmonic well or periodic box Calogero 

thermodynamics narrow down to a set of defining equations analogous to 

We will show that the thermodynamic of the Calogero model in a harmonic well can 
indeed be rewritten "a la TBA" in terms of a discretized function $ which will encode the 
statistical Calogero exclusion between different discrete harmonic energy levels and which 
will, as it should, reproduce, in the thermodynamic limit u — > 0, — k 2 ) in ([5]). Not 
surprisingly, the same conclusion will be reached for the LLL-anyon model in an harmonic 
well, whose thermodynamics will obey the same TBA equations as the Calogero-Moser 
thermodynamics. 

We will also argue that the same logic applies in the Calogero- Sutherland case, provided 
that a global shift of the bare quantum numbers is made in order to maintain a symmetric 
repartition of the dressed quantum numbers around zero. 
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Finally, we will look at possible applications of discrete TBA thermodynamics beyond 
the Calogero-Moser and harmonic LLL-anyon cases, by considering the Lieb-Liniger model 



in a harmonic well fH|, [19. This model is interesting because of its relevance to the 



description of one dimensional trapped Bose condensates [[[(J. We will show that the N- 
body groundstate energy is correctly reproduced at first order in perturbation theory by the 
discrete TBA equations, but corrections do appear at second order. 

II. THE CALOGERO CASE 

For a system with a discrete 1-body harmonic spectrum, the TBA equations (Q^) should 
rewrite quite generally as 

oo 

logZ = X>g[l + ze-*< n >], (7) 

n=0 

where the 1-body dressed energy e(n) should now be defined in terms of the 1-body Id 
harmonic spectrum (bare spectrum) e {n) = (n + 1/2)cj, < n as 

oo 

Pe(m) = /3e (m) - £ $ ni , n2 log[l + ze-^) (8) 

n2=0 

are just the discretized versions of (0,|). In (|), $ nii „ 2 has to be understood as acting 
on the free harmonic spectrum, i.e as acting on the power series in z obtained from (|8|) by 
expanding y{n) = 1 + ze' 13 ^ as 

oo 

V(n) = E VN{n)z N (9) 

N=0 

The lowest order terms of (|9|) are 

y (ni) = 1, 
y l (n 1 ) = e-^\ 

oo 

y 2 (ni) = £ e-^ n ^ nuna e-^\ 

n2=0 

OO , 1 s 

^ Y ni,n2 c ^ri2,ri3 ~ 2 n l> n 2 "1,713 y 



"2, ".3=0 
OO 



- o E e-^ (ni) $ nil n 2 e- 2 ^ (n2) , (10) 

2 n 2 =0 



where < ni, n 2 , • • • and the summation should be taken for all possible independant integers 
n 2 ,n 3 , ■ ■ -. 

If the TBA cluster coefficients obtained from expanding log Z = J2 b n z n 



oo 

b i = ^Vii n 

n=0 

oo 

6 2 = E 

n=0 

oo 

63 = E 

n=0 



2/2 (n) - (n) 



2/3(n) - yi(n)y 2 (n) + -yf(n) 



(11) 



have to match against the Calogero-Moser cluster coefficients &i 



e~' 3 " a -e-' 3 " __ 1 1 
(l_ e -2^)( 1 _ e -^) 2 l-e~ 2 ^ 



, . . ., $ni,n 2 should be defined as 



,n 2 -fnijnaV^) Pni,n 2 (<^ 1) 



(12) 



where P„ ljn2 (a) projects the two independant quantum numbers < ni, n 2 on dressed quan- 
tum numbers which, not surprisingly, obey exclusion statistics. More precisely, evaluating 
in ([7U^JT0|) expressions of the type 



J2 P ni,n 2 (oc)e 
n 2 =0 

P niiTl2 (a) amounts, n\ being given, to the shift 



-/3e(n 2 ) 



(13) 



n 2 ^ ni + n 2 + a, n 2 >0 



(14) 



and the summation over n 2 is replaced by the summation over n 2 > 0. In other words, in 
terms of the independant quantum numbers < n\,n 2 , denoting ni = ni,ri2 = n\ + h 2 , 
-Pni,n 2 (tt) means rt\ — > n 1; and n 2 — > n 2 + a, where < n x < n 2 are now bosonic quantum 
numbers. Therefore P ni)n2 (0) projects < ni,n 2 onto bosonic quantum numbers < ni < 
ri2, whereas P nijn2 (l) projects < rii,n 2 onto fermionic quantum numbers. Note that in 



(|T2|) substracting P ni:n2 (a = 1) is simply a matter of convention, i.e. as stressed above, a 
fermionic thermodynamical potential (|3|) with a spectrum which has to coincide with the bare 
spectrum when a = 1 -the Bose convention would yield $m,n 2 = Pni,n 2 (o ; ) ~ P ni ,n 2 ( a — 0). 
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More generally notice that ([12]) allows to rewrite (§) as 

y(m) - ^in^ Mn. + ^ + a) = 

1 ; nr 2 =o2/K+n 2 +i) 1 ; 

which can be viewed as the discretized version of (Q). 
Going one step further one gets 



DC 



II V( n i + ™2)= II y( n i + n 2 + l) + ze~ Pe ° {ni) n y(n x + h 2 + a) (16) 

n2=0 fi2=0 fi2=0 

which in turn, taken at n\ = 0, rewrites as 

oo oo 

Z = ze~^ J] Vfa + «) + II Vfa + 1) ( 17 ) 

h2=0 n,2=0 

Equation (|17|) can be interpreted by saying that either the first particle is in the ground- 
state at energy and then the next particle is in the energy level at least higher than 
(~ + a)u, or the groundstate is vacant and the first particle is in the energy level at least 
higher than (| + 1)c<j. One verifies recursively that Z in ( |TTD is identical to the grand partition 
of the Calogero-Moser model with iV-body spectrum 

N ( 1\ 
E N = ooY i ([n l + a(z-l)} + -J (18) 

where < n x < n 2 < . . .. In terms of the bare independent quantum numbers < ri\, n 2 , ■ ■ ., 
one has Hi — > rij + a(i — 1) = n^i + hi + a(i — 1) with hi > 0. This is indeed a BA "like" 
spectrum, i.e. in terms of the dressed quantum numbers = rij + a(i — 1), one has 
n',; = ni + aY^j^i — n'j). In the 2-body case, it indeed amounts to ri\ — > rii,n 2 — > 
rii + h 2 + a, i.e. to the the action of the projector P ni n2 (a) on the independent quantum 
numbers < ni,n 2 - 

Note that (0) implies that the iV-body partition function Zn obtained from (0) as 

oo 

Z N =l[y Nn (n), 0<N n , £ N n = N (19) 

n=0 n 

has, using (|T0|) to all orders, the simple factorized form 

oo 

Z N = PnunMPmM ■ ■ ■ P nN ^,n N ( a ) e -^(^+^(n 2 )+...+ eo (n N )) (20) 

ni,...,njv=0 



In particular in the 2-body case 

oo 

Z 2 -Z 2 \ Fermi = £ (P ni , n2 (a)-P ni , n2 (a = l))e-^^) + ^ 2 » (21) 

ni,ri2=0 

and thus in the thermodynamic limitQ uo — > 

^2 " Zaljrermi = (— f / dk.dk^kx - k 2 )e~^ + ^ (26) 

where $ is given in (f|). 

As far as the LLL-anyon model in a harmonic well is concerned, one finds that 
(0,|],[12D are unchanged, to the exception of the 1-body energy which now reads e (n) = 
(u t — uj c )n + uj c , where uj c = eB/2, u t = ^Juj 2 + u 2 , and the statistical anyonic parameter 
has to be understood as being —a, i.e. the screening regime where the flux <fi = —acf) (4> is 
the quantum of flux) carried by each anyon is antiparallel to the external magnetic field. 

One can easily convince oneself that in the thermodynamic limit u — > 0, both the 
Calogero and LLL-anyon TBA thermodynamics narrow down to 



2 By factorizing the center of mass, ( |2l|) rewrites as 

Z 2 -Z 2 \ Fermi = —^ J2 (e-^+i+a) _ e-M^i+D) (22) 

cm n I— — \ / 



Zblllll 2 n 2=2 i> 



and (H) as 

^2 - Z 2 | Fermi = (— ? / cOTe"^ / dk^(2k)e~^ (23) 
Since, in the thermodynamic limit 0] for the iV-th cluster coefficient, — ► —^=^/N, one infers 



that in the 2-body case 



2 sinh ^ 2tt 

and therefore one should have 

n 2 =2«>0 

a result that can be trivially checked by direct computation 



1 L C°° rk 2 

dKe' 13 — (24) 



J- ( e -M**+§+«) _ e -/Mn 2+ !+i)^ _^ L r dk<5>(2k)e-? k2 (25) 
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roo 

\ogZ= / de oPo {e ) ln(l + ze~^) (27) 
Jo 

where the dressed energy e(e G ) is implicitely defined a la TBA in terms of the bare energy 
Go as 

roo 

p e = p to - / de'$(e, e') ln(l + ze^ e ) (28) 
Jo 

and 

<t>(e,e') = (l-a)5(e-e') (29) 

Here, p (e ) is the 1-body density of states of the bare spectrum of the system considered, 
i.e. in the Calogero case the Id free density of states p (e ) = L/ (ny i /2e^), and, in the LLL- 
anyon case, the 2d LLL density of states, p (e ) = (BV/<f> )S(e — oo c ), where V is now the 
infinite surface of the 2d plane. 

In the case of the Calogero model in a periodic box -the Calogero-Sutherland model-, one 
can still propose (0) and but now the 1-body dressed energy e(n) should be defined in 
terms of the 1-body spectrum in a Id periodic box of length L, e G (n) = k 2 /2, with discretized 
momentum k = 2irn/L and n positive, null or negative integer. 

However, and contrary to the harmonic case, the very fact that the bare quantum num- 
bers in a periodic box can be of both signs lead to some adjustments. If one looks at the 
iV-body Calogero-Sutherland spectrum^ 

l / 2n, 9 J^f r , „ (iV-l)\ 2 . , 

E n = ^( T ) 2 E [\Pi + a(< " 1)] " a —^) ( 3 °) 

where ni < 112 < . . ., one finds that in terms of the bare quantum numbers rij — > + 
a(i - 1) - a(JV - l)/2 = + n { + a(i - 1) - a(N - l)/2 with hi > 0. This is quite 
similar to the Calogero-moser spectrum, up to a global shift, rij — > — a(N — l)/2, a 
iV-dependant periodic boundary condition adjustment insuring that the dressed spectrum 



3 with a BA spectrum u.[ = rii + ^ Sjyi sign(n'i — n'j) 
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remain symmetric around in order to minimize the iV-body energy. In the 2-body case 
U\ — > rii — a/2,n 2 — > n\ + h 2 + a/2, it amounts to the the action of the Calogero-Moser 
projector P ni>n2 (a) as given in ( |T2"D on the a priori two independent quantum numbers n\, n 2 , 
again up to the 2-body shift — > ni^— a/2. This being considered, altogether with the fact 
that the Calogero-Moser and Calogero- Sutherland models originate from the same model, 
up to a long distance regularisation, it is natural to take for both models the same TBA 
function (|l^) to obtain, in view of (|2~0l), the correct Calogero- Sutherland iV-body partition 
function, but in addition the shift n« — > rii — a(N — l)/2 has to be made a posteriori. 

At this point, one can remark that in all cases studied so far, the Calogero-Moser model, 
as well as the Calogero -Sutherland model up to periodic boundary conditions adjustments, 
and their thermodynamic limit, the Calogero model, the TBA functions $ nii „ 2 and $(£4 — k 2 ) 
are intimately related to the relative 2-boson density of states for the problem at hand. 
Indeed, in a harmonic well, the spectrum for a relative particle with bosonic statistics and 
interacting with a Calogero potential at the origin is 

e = C(J ( n+ I + a ) 0< n (31) 

with n even, i.e. with symmetric eigenstates under x — > —x. 

Let us first consider, in the thermodynamic limit, the Calogero model: when u — > 0, the 
relative 2-boson density of states reads 

Pa(e)-p a =i(e) = ^^5(e) (32) 

It rewrites in terms of the relative momentum k such that e = k 2 

1 - a 

P«(k) - Pa =i(k) = —5(k) (33) 

Now one has to map the relative 2-body momentum k on the "momentum" k 2 — k± the 
function $(k 2 — ki) is concerned with. Since k 2 — k\ = 2k, one gets for the density of states 
in terms of k 2 — k\ 

P^^ 1 ) - Pa=i(^y^) = (1 - - h) (34) 
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i.e. precisely (§) up to a factor 2n/L. 

When to is kept finite, the same logic applies: the relative spectrum fl3~T|) yields 



n — > n + a (35) 

One has yet to map the relative 2-body bosonic even quantum number n on the "quantum 
number" n 2 — rii that the function § ni ,n 2 is concerned with. For a given 2-body energy, i.e. 
for ri! + n 2 given, one has n 2 — n x = n -then the center of mass quantum number is 2ni, or 
n 2 — n x = n + 1 -then the center of mass quantum number is 2n x + 1, depending if n 2 — n x 
is even or odd. One finds that n — > n + a rewrites, in terms of the independent ni,n 2 as 
n\ — > ni, n 2 — > n\ +n 2 + a, where now n 2 = n or n 2 = n+ 1, i-e. any positive integer. Then 
( p5|) is indeed identical to the action of $„ b?l2 in (|12]). 

This is not a surprise, scattering 2-body phase shifts are known to be linked to the 2-body 



density of states via S-matrix arguments [18 



III. THE LIEB-LINIGER CASE 

In the Lieb-Liniger model in the thermodynamic limit, the same conclusion happens to 
be true. The model, defined as 

1 N d 2 

H N = ~^Y.JJ2 +C Y. 5 ^- X j) ( 36 ) 

i=l i i<j 

is solvable by Bethe ansatz fl4| and has a TBA thermodynamics |§ obtained from 



«fr- fa > = i (*» (37) 

It interpolates between the Bose (c = 0) and Fermi (c = oo) thermodynamics and describes 



particles with intermediate statistics ||15|| . For a relative particle with bosonic statistics 
interacting with a 5 potential at the origin the density of states is 

1 c 

p c (e) - p c=00 (e) = - — (38) 
which in terms of e = k 2 , k > rewrites as 
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1 c 

Pc(k) - p c =oc(k) = — — — (39) 

27T k 2 + \ 

Now, one has again to map k on the "momentum" k 2 — k\ the function $(&i, £2) is concerned 
with, i.e. hi — k\ = 2k, and since k 2 — fci can be either positive or negative, one gets for the 
density of states in terms of k 2 — k\ 

( k 2-k 1 ,k 2 -k x 1 2c 

Pc( ^— } " ^(-2—) = ^ (fr-fc^ + c* (40) 

i.e. nothing but fl3~T|), again up to a factor 2ir/L. 

If one follows the same line of reasoning which was operative in the Calogero-Moser case 
to obtain the discrete TBA function <E> niin2 (|12D from the relative 2-body spectrum ([H]), one 
might try, for the Lieb-Liniger model in a harmonic well, discrete TBA thermodynamics (0||) 
equations with a TBA function $m,n 2 deduced from the 2-body relative bosonic spectrum 



in a harmonic well [151. It rewrites as 



f _( n + i + l al ,ta„(^£||±|)) 0<» (41) 
with n even, i.e. as e = uj{n + ~ + f c (n)) withP] 

< fJn) = - arctan = { \ f 2 ' < 1 (43) 

vr \ v 2 v / 27Jr(^ + ^)y " 1 ; 

and interpolates between the relative 2-body bosonic (c = 0, fo(n) = 0, <7o( n ) — 1) an d 
fermionic (c = 00, /oo(n) = 1, g<x>(n) = 0) spectra. 

Therefore let us try for the Lieb and Liniger in a harmonic well the discrete TBA function 

^n lt n 2 = Pm,n 2 ( c ) ~ Pni,n 2 (c = Oo) (44) 

should be defined in terms of P ni ,n 2 ( c ) sucn that, n\ being left unchanged, 



^Equivalently, starting from the Fermi spectrum by rewriting f c (n) = 1 — g c (n) 

< o c (n) = - arctan — ^ < 1 (42) 

" V ^ TT ^ C r(^-2£M)y " V ^ 
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n 2 -> ni + n 2 + / c (n 2 ) (45) 

if n 2 > is even, and 

n 2 -> ni + n 2 + / c (n 2 - 1) (46) 



if n 2 is odd. Note again that substracting P ni n2 (c = oo) in fl44]) originates, as in the Calogero 
case, from the fermionic convention (obviously P ni) „ 2 (c = oo) = P m7l2 (a = 1).) 

It is easy to check that the 2-body partition function is reproduced by the discrete TBA 
equations 

oo 

Z 2 -Z 2 \ Fermi = J2 (P ni ,n 2 (c)-P nun2 (c = oo))e-^ n ^ n ^ (47) 

n i=n2=0 

and thus, in the thermodynamic limit uj —>■ 0, 

£ (e-flfia+i+Wfia)) _ ( c = ^ _, il /°° dk$(2k)e-P k2 (48) 

n 2 =2Z>0 27F 

where the function $ is given in fl37|), a result that can be checked by direct computation, 
order by order in 1/c. There are two independant dimensionless parameters, (3uj (thermo- 
dynamic limit) and ("coupling constant"). Clearly, for a given coupling constant y/]3c, 
looking at f c (n) = 1 — g c {n) in ( f4"3| , |42| ) , one has to consider, in the thermodynamic limit 



(3uj — > 0, the spectrum close to the Fermi point (c = oo), 

+ ... (49) 



from which (^) can be recovered, here at first order l/(^/J3c). 

Note also that in the thermodynamic limit for the relative spectrum, with {n + 1/2)uj — > 
k 2 , i.e. nu fixed, (|43D becomes 

2 c 

f c {k) = - arctan — (50) 

7T IK 



which is indeed reminiscent of the Lieb and Liniger BA spectrum [JT4 . 

What about the A^-body problem? A possible way to check the discrete TBA is to see 
if the perturbative TBA thermodynamics coincide with the exact (standard) Hamiltonian 
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perturbative thermodynamics |19j , which can be computed with the Lieb and Liniger Hamil- 
tonian from the Bose point c = (from the Fermi point c = oo standard perturbation theory 
is meaningless). Perturbation theory yields 

oo s+t fBuj 

\ogZ = \ogZ\ Bose + £ VP + . . . (51) 

V s % ,/sinh ^ sinh ^ sinh 



Let us now consider the y/]3c expansion from the TBA point of view, again from the Bose 
point. One has to consider (]43|) at first order in y^ftc 

m) 1 + ,52, 



vrv/2 v 7 ^ Z! 
and compute from the discrete TBA equations 

00 ze -/M n i+|) 00 x 
\ogZ = \ogZ\ Bose + £ , ^ $ ni ,n 2 ln- _ Mn2+1) |i + ... (53) 

ni=0 l — Ze pu,{ni + 2> n2= Q 1 — Ze P^Wi + i) 

where $ m n2 In — — — -tt|i means evaluating this expression at first order in sf^c using 

(|52|). One finds 



\ogZ = logZ\ Bose + j^-f^ (^ coth "I" + \/ coth T-J + ■ ■ ■ ( 54 ) 

which coincides with fl5T|) only in the limit /?a> — > oo, i.e. for a given cj, in the zero temper- 
ature limitf], i.e the groundstate. 

In fact, discrete TBA gives in the vanishing temperature limit direct information on the 
iV-body groundstate energy : in the Calogero-Moser case, it is obtained, in the bosonic 
based formulation, by restricting the discrete TBA equations 

OO 1 

\ogZ = log- —, — r, (55) 

m=0 1 ^ e 



5 The limit uj — ► oo for a given temperature is not considered here. In this limit all particles are 
confined at Xi = 0. But 5 interactions actually forbid this unless the effective coupling constant 
vanishes, which is precisely happening in the 2-body case (f43|). In other words the u — ► oo limit is 
the trivial bosonic limit. 
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/3e(m) = /3e (m) - £ $ ni>na log - _ 0i{n2) (56) 

ri2=0 

with $ nii n 2 = -Pni,n 2 ( Q; ) ~~ -Pni,n 2 ( Q; = 0) to the groundstate quantum numbers n\ — and 
n2 = 0. One obtains 

B - = w g + ^zi) a ) (57) 

In the Lieb and Liniger case a similar approach gives 

^ =w g + M /c(0) ) ,58) 

a result which is consistent with the groundstate energy at first order in c 

£ - = .g + fci) / (., ( )c+ ...) ,59) 

where /«(0) = 1^/2 stands for the first order term in the expansion of / c (0) in power 
of c. However, second order standard perturbation theory gives corrections to the Calogero- 
Moser like energy uN(N — l)/ c (0)/2. For example in the = 3 case 

Ef S - = u, + 3/« (0)c + 3 (/i 2 ) (0) - -L log c 2 + . . .) (60) 

where /j 2 ^(0) = — stands for the second order term in the expansion of / c (0). 



IV. CONCLUSION 

We have shown how the Caloger-Moser thermodynamics can be rewritten in terms of 
discrete TBA equations. In the Calogero-Sutherland model, the same TBA equations were 
shown to be operative, up to a global shift of the bosonic quantum numbers. Since the 
Lieb-Liniger model shares common features with the Calogero model -BA solvability, TBA 
thermodynamics in the thermodynamic limit, intermediate statistics- it might also have, 
when considered in a harmonic well, a discrete TBA thermodynamics. We tried to illustrate 
this point of view by proposing discrete TBA equations for the harmonic Lieb and Liniger 
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model in analogy with the Calogero-Moser TBA thermodynamics. However the groundstate 
energy shows deviations from this TBA framework at second order in perturbation theory. 

We leave to a further study to find analytical or numerical ways to improve and give 
a stronger basis to the discrete TBA thermodynamics for the Lieb-Liniger model, and in 
particular extract a useful information on the groundstate for a given density of particles. 
It would however certainly be interesting to understand more in detail the zero temperature 
limit of a system which is supposed to describe the physics of Id Bose Einstein condensates 
in harmonic traps. 

Acknowledgements: One of us (S.O.) would like to thank S. Matvenkoo for useful dis- 
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